Forbidden island heights in stress-driven coherent Stranski-Krastanov growth 
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The observed height distribution of clusters obtained in strained epitaxy has been often interpreted 
in terms of electronic effects. We show that some aspects can be explained classically by the interplay 
of strain and edge energies. We find that soft materials can transform directly from monolayer 
into thicker islands by two-dimensional (2D) multilayer nucleation and growth. There is a critical 
thickness decreasing with the force constant. Thinner islands are thermodynamically forbidden, due 
to the insufficient stress relaxation upon clustering particularly under tensile stress. At sufficiently 
large misfits the barrier for 2D multilayer nucleation is significantly smaller than the barrier for 
subsequent single-layer nucleation. The effects are found to be quantitatively reasonable and offer 
a plausible explanation for the absence of thin islands and 2D growth of ffattop islands usually 
attributed to quantum size effects. 



PACS numbers: 68.35.Md, 68.43.Hn, 68.55.Ac, 68.65. Hb 
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Nanostructures are very promising for optoelectronic 
and magnetic applications. For efficient operation, the 
shape, size and thickness distribution of small clusters 
are important parameters. Therefore it is crucial to un- 
derstand the factors which control them. The epitaxy 
of metals on semiconductor surfaces at low temperatures 
(130K - 180K) has been intensively studied in the last 
years 0, [3, y, H- Some important observations are: 
(i) flat-top Pb islands with steep edges and a preferred 
height of 7 monolayers (ML) grow on the wetting layer 
on Si(lll)7 X 7 [l|, "11, i, H, (ii) islands with thicknesses 
from 1 to 3 MLs are never observed 0, H, Q ; (iii) fiat-top 
Pb islands are preceded by pyramidal or dome-like clus- 
ters; clusters thinner than 3 MLs are never registered Q; 

(iv) 2-ML thick flat-top Ag islands on Si(lll) increase 
linearly in size with increasing coverage whereas ML is- 
lands preserve a nearly constant size of about 5 nm^ Q ; 

(v) flat-top islands with a preferred height grow later- 
ally without thickening 

[II, SH; (vi) Vertical growth of 
Pb/Si(lll) takes place by bilayer increments [tI. [loj. 

The above observations were explained in terms of the 
energy lowering due to electron confinement and spilling 
of charge through the metal-semiconductor interface by 
Zhang, Niu and Shih, who coined for this reason the term 
"electronic growth" However, classical effects asso- 
ciated to strain relaxation are expected to contribute, as 
assumed to explain the fact that the increase of the aspect 
ratio of flat-top Pb islands is mediated by the formation 
and growth of strips (or rings) around the outer edges of 
the islands The aim of this paper is to show that 

some of the observations listed above, in particular (ii), 
(iii), (iv) and (v), can be explained classically in terms 
of the interplay of strain and edge energies. 

Two mechanisms have been invoked to address the in- 
stability of planar growth against clustering. The first is 



the nucleationless development of instabilities of a certain 
wavelength that evolve into faceted three-dimensional 
(3D) islands [l^l- The second is a nucleation behavior 
due to the competition of surface energy and strain re- 
laxation [l3| . Studies of the total energy per atom of is- 
lands of different heights as a function of their total num- 
ber of atoms have shown that intermediate states with 
thicknesses increasing by ML steps are stable in separate 



consecutive intervals of volume |15l . Il6l | . The rearrange- 
ment of mono- to bilayer islands was found to be a true 
nucleation process in compressed overlayers, in the sense 
that a small critical nucleus of the second layer is ini- 
tially formed and then grows further up to the completion 
of the transformation. The associated energy displays a 
maximum at some number of atoms in the second level, 
which becomes very small for high enough values of the 
lattice misfit [l3|- However, bilayer islands of expanded 
overlayers tend to require unrealistically high values of 
the lattice strain and of the island size to become stable 
against ML islands. Their transformation curves show 
an increase of the energy up to large second layer cluster 
sizes, making nucleation unlikely Also compressed 

overlayers of softer materials, or at lower misfits, are not 
expected to show nucleation behavior through this mech- 
anism. We therefore set out to study the mono- to mul- 
tilayer island transformation by using a simple atomistic 
model in 2 -f 1 dimensions. 

We consider fcc(IOO) crystallites with the shape of 
truncated square pyramids on a rigid layer of the same 
material with lattice constant a. In our model, atoms 
interact through a pair potential 



V{r) = 



II — V 
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that in spite of its simplicity includes all necessary fea- 
tures to describe real materials (its strenght and anhar- 
monicity are governed by the constants ^ and v [lit). 
The equilibrium atom separation is 6, so that the lattice 
misfit is given by e = (5— a)/a. We consider interactions 
only in the first coordination sphere. For the study of 
mono- to multilayer transformations, we assume model 
processes in which atoms are detached one by one from 
the edges of an initial square ML island, placed on top 
of it at the center and arranged in a (X — 1)-ML thick 
structure as compact as possible up to the formation of 
a X-ML truncated pyramid. For every intermediate con- 
figuration, the cluster is allowed to relax by varying iter- 
atively the atomic positions until all the forces fall below 
some negligible cutoff value. The energy change associ- 
ated with the transformation at a particular stage is given 
by the difference between the energy of the incomplete 
multilayer island and that of the initial ML island. 
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FIG. 1: Energy per atom of mono- and multilayer islands 
as a function of the total number of atoms. The misfit e 
amounts to -7% and fi — 21/ — 26. The insert shows the 
dependence of the critical size Ni-x (beyond which 1-ML 
islands become unstable against X-ML islands) on the force 
constant 7. The threshold thickness X is denoted by the 
numbers at each point. 

Figure [T] shows the total energy per atom of expanded 
islands of different thicknesses as a function of the to- 
tal number of atoms. It is seen that monolayer islands 
remain always stable against bilayer ones. At a critical 
size trilayer islands become stable and, with sizes 

increasing further, the ground state shifts to islands of 
thicknesses increasing in ML-steps. The insert of Fig. [1] 
illustrates the dependence of the critical size on 
the force constant 7 {^ivVo in our model) for a fixed 
value of the misfit (e=-0.07). The threshold thickness be- 
low which multilayer islands are energetically unfavored 
steeply increases together with the critical volume with 
increasing softness (decreasing 7) of the epitaxial over- 
layer. 

The stability regions of mono- and multilayer islands 
for both compressed and expanded overlayers are visu- 
alized in the phase diagram of Fig. [21 as a function of 
the number of atoms in the islands and the bulk strain 
energy per bond £ = 0.5je^a^, which turns out to be the 
relevant scaling factor. As seen, ML islands are favoured 
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FIG. 2: Phase diagram showing the stability ranges of islands 
of difi^erent heights (given by the numbers in ML's) in coordi- 
nates of the total number of atoms and the bulk strain energy 
per bond £ = O.Sje^a^, for positive and negative values of the 
misfit and for fi = 2u potentials. 



at small number of atoms. For decreasing strain ener- 
gies, the transition takes place for larger island sizes and 
to thicker islands. Expanded overlayers require larger 
strain energy to become unstable against clustering. Is- 
lands thinner than a certain number of layers are for- 
bidden for energetic reasons. Result are shown for the 
particular class of potentials with /i = 21/; other choices 
of fj, and 1/ give points very close to those displayed. Dif- 
ferent potentials (e.g. Lennard-Jones) are expected to 
give extremely similar results. 

To compare with experiments, we estimate the strain 
energy per bond £ from the force constant 7 of the over- 
layer material (in turn calculated as 7 = Eh/2{\ — i/p), 
with E the Young modulus and i/p the Poisson ratio). 
Expressing £ in units of the bond energy Vq of the cor- 
responding material, we get 0.15, 3.08 and 2.27 for the 
cases of Pb, Ag and Al on Si(lll), or 0.05, 1.5 and 1.3 
eV, respectively. Thus, a transition for Pb/Si from 1 ML 
to 8 or 9 ML for islands of several thousands of atoms 
is expected, in reasonable agreement with experiments 
showing a preferred height of 7 ML and scarce presence 
of thinner islands [l|, i, S, H, S 13, S- For Ag and Al, 
the unrealistically high values of £ (due to the failure 
of the harmonic approximation for such high values of 
the negative misfit) predict no absence of thin islands, 
as observed [1, [l^ (see below). These estimations con- 
sider (100) coherent islands while most experiments are 
done on Si(lll) and the islands are probably relaxed. 
However, they may retain approximate validity because 
the strain energy relaxation is partially balanced by the 
cost of disregistry. Furthermore, calculating similarly for 
the systems In and Cd on Si (to our knowledge not yet 
studied), we obtain 0.17 and 2.65 Vq, respectively, thus 
predicting for In a similar behaviour as Pb (thin islands 
forbidden) in contrast to Cd. 

The results above point to a possible mechanism of 
2D-3D transformation different from the consecutive for- 
mation of bilayer, trilayer, etc. islands i7,23|: the direct 



nucleation of multilayer clusters whose thickness depends 
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on the values of the force constant and the lattice misfit. 
We have studied this mechanism: Fig. [3] shows transfor- 
mation curves from 1 to 3 ML islands for two negative 
values of the misfit, -0.07 and -0.12. Atoms are detached 
from the edges of the initial monolayer island and incor- 
porated to the double steps of the bilayer island growing 
on top. The low-misfit curve is similar to the layer-by- 
layer curve shown in Fig. 5b of Ref . J/7|. The energy tends 
to increase all the way (the departure from monotonic be- 
havior depends on the exact atomistic processes related 
to filling and depleting of atomic rows during the trans- 
formation process) and shows at the very end a sudden 
collapse due to the disappearance of the single and dou- 
ble steps to produce low-energy facets. On the contrary, 
the large-misfit curve shows in the beginning a nucleation 
behavior with a bilayer nucleus consisting of 22 atoms. 
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FIG. 3: Transformation curves representing the energy change 
in units of the bond energy Vo as a function of the number of 
atoms transferred to the upper level for (a) e = —0.07, and 
(b) e — —0.12. The number of atoms in the initial monolayer 
island (365 ~ 19 x 19) gives a complete truncated trilayer 
pyramid (12 x 12 + 11 x 11 + 10 x 10) and n = 2v = 26. 

We consider first transformations from mono- to bi- 
layer islands. Two processes accompany the island thick- 
ening: first, the distance I between the upper and lower 
steps decreases. Since the step repulsion energy increases 
only as [2lj, it is neghgible except very close to 
the end of the tranformation. The second is the es- 
sential process, the increase of the total length of the 
steps, which in our model of square islands is given by 

^iVo), with iVo and N2 being 



AL = 4(ViVo-iV2-fV^- 
the number of atoms in the initial ML island and in the 
incomplete upper layer, repectively (note that, as a first 
approximation, there is no change in surface or interface 
energies, since there is no variation of the out-of-plane 
coordination of the atoms, and islands and wetting layer 
are composed of the same material). This process has 
three energetic contributions associated: the first is an 
increase in step energy, understood as the energetic cost 
of the unsaturated bonds, due to the reduced coordina- 
tion at steps. The second is the relaxation of strain at 
steps and the third is the reduction of interlayer adhe- 
sion at steps due to the edge atoms climbing up in the 
potential valley created by their neighbors underneath. 
Only the strain relaxation favors the process of cluster- 



ing llZl, [24! . The other two contributions oppose it, in 
addition to the step repulsion energy. 

The sign of the misfit strongly affects the strain re- 
laxation at steps. For both signs and for large enough 
islands, the bonds are almost completely relaxed at the 
edges and strained at the center. However, for interme- 
diate values of island's size and misfit (3 - 7%) [H,!!!], 
the region with partially relaxed bonds extends consid- 
erably further in the case of compressed islands. Thus, 
as a result of anharmonicity, strain relaxation at steps 
is stronger for positive than for negative misfits of the 
same absolute value. The reduced coordination at steps, 
i.e. the number of unsaturated bonds, The number of 
edge atoms obviously does not depend on the sign of the 
misfit, while the reduced adhesion to out-of-plane neigh- 
bours is expected to depend weakly, because the inter- 
action of an edge atom with its neighbours in the layer 
underneath necessarily involves attractive and repulsive 
forces in both cases. In compressed islands the strain re- 
laxation at steps overcompensates the increase of the step 
energy (the total step length increases very fast in the 
beginning of the transformation). The balance turns for 
an island size which becomes very small for large enough 
values of the misfit. This results in nucleation-like behav- 
ior [itI . In expanded islands the weaker strain relaxation 
at steps requires larger values of the misfit to overcome 
the remaining contributions. 

We consider next the instability of mono- against mul- 
tilayer islands at small values of the force constant of the 
material. Owing to the weak strain relaxation at steps, 
the strain energy per atom decreases very weakly with 
thickness. A useful picture is to consider n-ML islands 
in coherent epitaxy as 1-ML islands but with n times 
stronger bonds, i.e. n times "stiff er" [2^. If the force 
constant of the material is small, the gain in strain en- 
ergy upon the formation of ML islands on the first ML 
will be also small. Instead, bilayer islands will become 
energetically favoured with increasing number of atoms, 
i.e., the ground state will evolve from mono- to trilayer 
islands. These, being effectively "stiffer" , relax the strain 
energy at steps more efficiently, overcoming the energetic 
cost of the bilayer step. For smaller values of the force 
constant tetralayers and even thicker islands will become 
stable against monolayer islands when the total size in- 
creases. This behavior is found also in compressed is- 
lands, shifted to smaller values of the force constant due 
to the anharmonicity of the interatomic potential. 

For both compressed and expanded overlayers, the de- 
crease of the strain energy per bond leads to an increase 
of the critical thickness X at which the corresponding 
island becomes energetically stable against a ML island 
of the same size. A nucleation-like behavior is expected 
for this 1-X ML transformation. This basic result of our 
study, summarised in Fig. [21 implies that islands with 
insufficient thickness, i.e. with insufficient relaxation of 
strain energy at steps, will be forbidden to form for ther- 
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modynamic reasons. This can explain the absence of Pb 
islands thinner than 4 ML [1,0, Sll]. 
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rectly into thick islands by 2D multilayer nucleation 
and growth. Islands thinner than the critical thickness 
are thermodynamically forbidden due to the insufficient 
strain relaxation at steps, particularly when under tensile 
stress. At sufficiently large misfits the barrier for 2D mul- 
tilayer nucleation is significantly smaller than the barrier 
for the subsequent single-layer nucleation. The above ef- 
fects due to the interplay of step energy and strain relax- 
ation are found to be quantitatively reasonable and offer 
a plausible explanation for the absence of thin islands 
and the 2D propagation of flat-top islands in epitaxial 
growth of metals on semiconductors. 



FIG. 4: Nucleation barriers in units of Vb as a function of 
the absolute value of the negative misfit, for mono- to trilayer 
islands. The number of atoms in the critical nucleus is given at 
each point. Also shown is the barrier for monolayer nucleation 
on top of the trilayer island. The total number of atoms is 
365, and ^ = 2v = 26. 
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Figure [5] shows a further important result. The mono- 
tonic decrease of the nucleation barrier height with misfit 
for the 1-3 ML transformation is characteristic of nucle- 
ation behaviour fv^. For small misfits, the barrier for 
the 3-4 ML transformation is smaller than for 1-3 ML, 
implying that, if 1-3 is possible, so will be also 3-4 or even 
higher. However, at larger absolute values of the misfit 
the barrier for formation of a bilayer nucleus on the first 
ML is significantly smaller than the barrier for forma- 
tion of a ML nucleus on the trilayer, so that the growth 
of the 4th atomic level might be kinetically inhibited in 
this case. A similar intersection of the curves (not shown) 
is found for positive misfits, shifted to smaller absolute 
values. Thus, for large absolute values of the misfit the 
strain relaxation at the ML cluster growing on the trilayer 
island is not strong enough, overcoming the increase in 
step energy, to result in a small barrier height (this is 
essentially the same behavior as in the 1-2 ML transfor- 
mation discussed above) Therefore, at sufficiently large 
misfits the growth in height can be strongly inhibited due 
to kinetics. For example, a factor 2 in the barrier height, 
as given by FigU] for high misfits reduces the nucleation 
rate several orders of magnitude at room temperature 
or below and 2D lateral growth will be preferred, as in 
fact observed 0, 0, Q ■ In addition, the side walls act as 
better traps for adatoms than the flat top surface. The 
lateral growth by rearrangement of mono- to multilayer 
islands is strongly supported by the observation that 2- 
ML thick flat-top Ag islands on Si(lll) increase linearly 
in size with increasing coverage whereas ML islands pre- 
serve a nearly constant size of about 500 Q, which 
can be interpreted as the critical size iVi_2- A clear evi- 
dence for a transformation process is the rearrangement 
of 2 to 3-ML Fe islands on Cu3Au(001) upon annealing 



at 400K [2i,[27|. 



In conclusion, we find that materials with weak inter- 
atomic bonding will transform from ML islands di- 
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